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Abstract 

Charmless B decay modes B —> 7T7t, ttK and KK are systematically investigated with and 
without flavor SU(3) symmetry. Independent analyses on 7T7 t and irK modes both favor a large ratio 
between color-suppressed tree (C) and tree (T) diagram, which suggests that they are more likely 
to originate from long distance effects. The sizes of QCD penguin diagrams extracted individually 
from 7 T 7 r, nK and KK modes are found to follow a pattern of SU(3) breaking in agreement with the 
naive factorization estimates. Global fits to these modes are done under various scenarios of SU(3) 
relations. The results show good determinations of weak phase 7 in consistency with the Standard 
Model (SM), but a large electro-weak penguin (Pew) relative to T + C with a large relative strong 
phase are favored, which requires an big enhancement of color suppressed electro-weak penguin 
(Pew) compatible in size but destructively interfering with Pew within the SM, or implies new 
physics. Possibility of sizable contributions from nonfactorizable diagrams such as IT-exchange 
( E ), annihilation (A) and penguin-annihilation diagrams(PA) are investigated. The implications 
to the branching ratios and CP violations in KK modes are discussed. 
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I. INTRODUCTION 


In the recent years, the two A-factories have succeed in steadily improving the measure¬ 
ments of hadronic charmless B decays. At present, all the branching ratios of B —> 7 T 7 t and 
7iK modes have been measured with good accuracy. The large direct CP violations have also 
been established in 7 r + 7 r - and tt + K~ modes jl]]. Their implications have been reported in a 
recent short paperIt has been shown that the weak phase 7 can well be determined to 
be consistent with the standard model, it prefers a relative large electroweak penguin with 
a large strong phase and also favors an enhanced color-suppressed tree diagram. In this 
longer paper, we shall provide a more detailed analysis with including subleading diagrams 
and their implications to KK modes as well as paying attention to SU(3) broken effects. 


It is of interest to note that the signs of the direct CP violations, if confirmed bvthe 
future experiments would agree with the results from perturbative QCD approach 3, dj 
while posse a challenge to the naive factorization (b| and QCD factorization calculations 


__ j These impressive new results have triggered a large amount of theoretical efforts 

to understand the strong interaction dynamics of those decays 00,131110 , extract 
the weak phases in the Cabbibo-Kobayashi-Maskawa (CKM) matrix 00 and explore 
new physics 00 00 , 000 . 

Making use of the flavor topology of the decay amplitudes and the approximate flavor 
SU(3) symmetry, one can describe those decay modes in terms of several independent quark 
flavor flow diagrams, such as tree diagram (T), color-suppressed tree diagram (C), QCD 
penguin diagram (V), electroweak penguin diagram ( Vew ) an d color suppressed electroweak 
penguin diagram ( Vew ) e ^ c - dt then follows from the hierarchies of the Wilson coefficients 
and the CKM matrix elements that the B —» mr modes are T dominant while B —> ttK 
modes are V dominant. Since C is color suppressed, one expects that the hierarchical 
structures among those decays should be 


2Br('K°TT °) -C Br(7T + 7T ) « 2Br(-K 7 T°), 


( 1 ) 


and 


Br(TT + K ~) ~ Br(7T~K °) ~ 2Br(7T°K °) ~ 25r(7r 0 A"”), (2) 


respectively. 

Note that the above relations follow from a purely short distance diagrammatic de¬ 
script ion which could be misleading in the presence of large final state interactions (FSIs) 
(25l l26l 1271. At present, they are not favored by the experiments. The current world aver¬ 
age data (lL l28lj listed in TablH show big enhancements of i?r( 7 r° 7 r°) and Ar( 7 r _ 7 r°) relative 
to Br( K 7i + 7i~) and suppression of Br(7T + K~) relative to 2Br{ir°K Q ) and Br(n~K°). The 
numerical values of these relative ratios are given by H3 




Roo 


2 Br(7T~n°) 
Br(7T + 7T~) 
^ Br(n°7i 0 ) 
Br( 7 T+ 7 T - ) 


• — = 2.2 ±0.31, 
t b + 

= 0.67 ±0.14, 


(3) 
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Br(x 10 6 ) 

acp 

5 

7r + 7r 

4.6 ± 0.4 

0.37 ± 0.11 

-0.61 ± 0.14 

O 

O 

1.51 ±0.28 

0.28 ±0.39 


O 

1 

5.5 ±0.6 

-0.02 ±0.07 


n + K 

18.2 ±0.8 

-0.11 ±0.02 


n°K°(K s ) 

11.5 ± 1.0 

-0.09 ±0.14 

(±0.34 ± 0.28) 

1 

0 

24.1 ± 1.3 

-0.02 ± 0.034 


1 

O 

12.1 ±0.8 

0.04 ± 0.04 


K+K~ 

K°K° 

1 in+0.42 

1 - ly -0.37 



K~K° 

< 2.4(1.4518;^) 




TABLE I: The latest world averaged data of Charmless B decays 0, ;2£]. 


and also 


Rn 

R 

R 2 


Br(n + K ) 

2Br(n°K°) 

Br(n + K~) 

Br(n~K°) 

Br(n + K ~) 

2Br(n°K~) 


= 0.79 ±0.08, 


t b + 

r B o 

t b + 

Tb° 


= 0.82 ±0.06, 
= 0.81 ±0.06. 


( 4 ) 


The above five ratios characterize the puzzling patterns of the latest data and may provide 
insights on the strong dynamics of heavy quark decays or possible new physics beyond the 
Standard Model (SM). 

The large value of i?oo forces the C to be large, which is a challenge to theory. Various 
ways to explain large i?oo with reasonable values of C/T involve an enhanced IV-exchange 
diagram (£1), a large QCD penguin contribution corresponds to w-quark loop or large final 
state interactions (FSIs) which involves DD^ intermediate states and quasi-elastic mixing 


between n + n~ and 7r°7r° modes 13] • The recent SCET calculations also supported a large 
C/T J2f|. Note that the tht and nK modes differ in flavor topological structure while FSIs 
are flavor blind, the two kind of effects can in principle be distinguished by a separated 
study of these two sets. FSI will lead to large C in all decays modes. Furthermore, it should 
enhance Pew relative to Tew in a similar manner. 

In the nK modes, it is well known that the suppression of R n is more relevant to the 
electro-weak penguin sector, as in nK modes T and C are greatly suppressed by small CKM 
matrix elements, In the SM, from the isospin structure of the effective Hamiltonian, the 
ratio between electroweak penguin and tree diagrams are fixed through SHE HQ 


tdSM _ 

rr ew — 


EW 


+ P& 


EW 


3 C 9 ± C' 


10 


T+C 


2 Ci ± C 2 


= (1.35 ±0.12) x 10 


-2 


( 5 ) 


for nn modes. Where T, C, Pew an d Pew are diagrams with CKM matrix elements 
factorized out which will be discussed in detail bellow. C*s stand for the short distance 
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Wilson coefficients at the scale of /i ~ mb- This relation is free from hadronic uncertainties 
and survives under elastic FSIs and inelastic FSIs through low isospin states such as B —> 
DD S —> 7T7t(A"). It also predicts the direct CP violation in 7r _ 7r° modes to be vanishing. 
Using flavor SU(3) symmetry, this relation also holds for nK modes. Thus it can directly 
confront the experiments and allows us to explore the new physics in hadronic charmless 
B decays. It is of interest that the charmless B decay data indeed imply the violation of 
Eq. ®. The possibility of larger isospin / = 2(3/2) amplitudes violatin g Eq . ltHl) in iminK) 


and recently discussed in 


mtudes violatin g Px 
Refs. |H, H E E,' 



modes was found in Ref.|3£ 
with updated data. In a recent analysis, an enhancement of a factor of two was obtained 
through a direct global y 2 analysis using flavor SU(3) symmetry^. 

Although it is too early to draw any robust conclusion, it motivates us to take a closer look 
at the electroweak penguin sector within and beyond the SM. Note that in these analyses 
on large Pew, the effects of Pp jW are often assumed to be small, which is conceptually not 
appropriate as P%w * s directly involved in Eq. ©• Furthermore, it provides a cancellation 
to the low isospin / = 0(1/2) part of Pew- Thus its effects could be significant. 

The suppression of R may require significant contributions from subleading diagrams 
such as annihilation diagram A or color-suppressed electro-weak penguin diagram Vew- 
Considering the fact that A contributes to n~K° and n 0 K~ in the same way, namely they 
have the same A — V interference, one expects that an enhancement of A with appropriate 
strong phase can suppress simultaneously R and R 2 while their effects would cancel in some 
extent in their ratio. The current data give 


Rh = 


2Br(n 0 K~) 
Br(ir-K °) 


= 1.0 ±0.08, 


( 6 ) 


which agrees well with this conjecture. However an enhancement of A will lead to significant 
consequences to KK modes, especially K~K °, as it is not suppressed by CKM matrix 
element like in the ttK modes. It is expected that a strong constraint on A will be found 
once this decay mode is experimentally observed. 

There already exists a number of global \ 2 analyses on B —> 7nr,iiK and KK systems 
based on flavor SU(3) symmetry Il2l . lOul . |4(1 14ll . j42j. I±1 l44l | . But to trace back the origins of 
the above mentioned 7T7 t and nK puzzles, separate y 2 analysis are urgently needed. Further¬ 
more, the SU(3) breaking scheme dependences are not carefully examined in the previous 
analyses, which may lead to different results in the literature. Finally, the contributions 
from subleading diagrams such as P^w , £■> Fl and penguin induced annihilation diagram V A 
which could play important roles in understanding the current data are often neglected. 

The purpose of the present paper is to make an up to the date investigation on charmless 
B decays, following a strategy that first applying y 2 analysis on 7T7 t, nK and KK modes 
separately, then connecting them through flavor SU(3) symmetry and discuss the SU(3) 
breaking scheme dependency. After obtaining reasonable values of the dominant amplitudes, 
we then discuss their implications to A" Jl modes with sub leading diagrams such as Vew, £, -A 
Va etc. 

Our main observations are 


• Independent fits on mr and ttK modes without SU(3) symmetry both favor a large 
ratio between color-suppressed tree (C)and tree (T) diagram, which disfavors the ex¬ 
planation of large nonfactorizable VF-exchange diagrams (£). The extracted QCD 
penguin diagrams from 7T7 t, nK and KK show a clear signal of SU(3) breaking and 
the breaking pattern is in agreement with naive factorization. 
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• Global fits for 7 T 7 t, nK and KK modes show good determinations of weak phase 7 in 
agreement with the standard model and prefer a larger electro-weak penguin ( Vew) 
relative to T + C with a large strong phase when P BW is neglected. The results are 
found stable among various SU(3) breaking schemes. The current data favor a SU(3) 
breaking scheme in which all the amplitudes for ttK are greater by a factor of fx/f-n 
motivated from factorization. 

• An enhancement of Pew with destructive interference to Pew provides a alternative 
explanation to the small R n within the SM. The suppression of R can be partially 
explained by an enhanced annihilation diagram A. The Va provides a source of SU (3) 
breaking in strong phases. 

• The subleading diagrams may lead to significant CP violations in KK modes. For 
typical value of A and V, the direct CP violation in K~K° can reach ~ 0.4. 

This paper is organized as follows. In section El the basic formulas for diagrammatic 
decomposition are presented. In 11111 we extract the parameters such as weak phase 7 and 
the decay amplitudes from 7 T 7 t, 7 xK and KK modes separately. In section ITV1 we combine 
them in three different scenarios of SU(3) symmetry. One is that all the amplitudes in nK 
modes are rescaled by a factor of fx/fn motivated from the native factorization. An other 
one is that only the tree diagrams are rescaled by this factor while the rest of the amplitudes 
remain the same in SU(3) limit. The last one is the strict SU(3) limit. In section 0 We 
consider the contributions from various subleading diagrams and extract their typical values. 
In section eh The implications to the KK modes are discussed. The possibility of finding 
large direct CP violations is indicated. We conclude in section IV 111 


II. DIAGRAMMATIC DESCRIPTION 


We use the following definitions for branching ratios and direct CP violations 


Br 

acp 


^(l^| 2 +l^| 2 ), 

I^| 2 -I ^| 2 

l^ + l^l 2 ’ 


(7) 


where A(A) stands for B°(B) or B + (B~) decay amplitude, r is a phase space factor, r = 1 
for neutral final states and r = t b +/tbo = 1.086 ±0.017 for charged final states. The mixing 
induced CP violation parameters S and C are introduced through the time-dependent decay 
rate difference 

_ r(fi° -> f C p) - r(a° -> fcp) 
aCP{ ’ r(B« f CP ) + r(B° fcp) 

= S ■ sin (A ms ■ t) — C ■ cos(A m B ■ t), (8) 


with fcp denoting final states with definite CP parities. A m B is the neutral B meson mass 
difference. The two parameters can be written as 


2Im\ _ 1 — |A | 2 
1 ± |A | 2 ’ “ 1 ± |A | 2 


—O-CPi 


(9) 
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with 


A = e~ 2 “>^ (10) 

in the SM. 

Using the phase definitions of B~ = (—ub),B° = ( db ), K~ = (— us),K° = ( ds ) and 
7r + = (ud), 7T° = (dd — uu)/ 2,7r~ = (— ud ), one arrives at the following diagrammatic 
decompositions for 7T7 t modes |45|, |46|, [47[ 


^4(7T + 7T ) — — (T + £ + V + Va + 7^Pew)i 

^(7r°7T 0 ) = ——^=(C — £ — V — Va + ’Uew + V EW )i 

A(A~) = --h(T + C + ~Pew + V C EW ). (11) 

Similarly, the ttK modes are given by 

A(* + K~) = -<r + v' + ^V% w ), 

,4(V J /\ ) — - V + V' EW + -V EW ), 

= A + v- l -v c EW , 

A(t°K-) = ~(T + C' + A' + V + ?V + |pg w ). (12) 

The amplitudes for nK modes are marked by a prime, which equal to the unprimed ones 
for 7T7T modes in flavor SU(3) limit. The KK modes are given by 

A(K+K~) = -(£" + V[), 

A(K°K°) = V" - + V" A , 

A(K~K°) = V" - l -V c EW + A", (13) 


where the subleading diagrams such as color-suppressed electro-weak penguin (V EW ), W— 
exchange diagram (£), annihilation diagram A and penguin-induced annihilation diagram 
(Va) are included. In the above formulas, the penguin exchange diagram ( Ve ) are ab¬ 
sorbed into penguin diagrams and the electroweak and color suppressed electroweak penguin 
exchange diagrams are neglected. 

We start with independent analyzes on 7T7 t, ttK and KK modes. In the first step, all 
subleading diagrams such as V EW , £, Va and A are switched off for the reason of simplicity. 
Their effects will be investigated in detail in section IIVI and El To consistently include 
the experimental errors of the data, the y 2 method is adopted for extracting the decay 
amplitudes. The definition of y 2 reads 


= S 

i 


r""(<\i)i - fi 


exp 


Vi 



(14) 
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where fl heo are the theoretical values of observables fi(i = 1 ,m) and a,j(j — l,n) are the 
to-be-determined parameters. f^ xp and cr i are the experimental central values and errors. 
The best-fit of the parameters correspond to the minimum of the y 2 function which satisfies 
a x 2 distribution with degree-of-freedom(d.o.f) m — n. For the experimental data we take the 
values listed in Tab|I]which are the weighted average of CLEO, Babar and Belle collaboration 
result sl28j |. Other major parameters used in the fits involve the CKM matrix element of V ub 
pH. !49l]and U c &|50j. In the numerical calculations we take the following values [28] 

V cb = 0.04 ± 0.02, V ub = (3.9 ± 0.68) x 10 -3 , (15) 

and the SM value of (Hi , [52.j 

sin 2/3 = 0.73 ±0.037. (16) 

All the Brs are written in units of 1CU 6 , and the angles are in gradient and arranged in the 
range (— n, ±7r). 


III. ANALYSIS WITHOUT FLAVOR SU(3) SYMMETRY 

A. 7T7T modes 


The hierarchies in the decay amplitudes are controlled by both the Wilson coefficients 
and the CKM matrix elements. As the sizes of the CKM matrix elements are better known, 
it is helpful to factorize them out so that all the hadronic amplitudes in 7T7 t, nK and KK 
etc have the same hierarchical structure. Thus we shall use the following parameterizations 
for 7T7r modes 


A{ 7r + 7r ) = — 


K(T + E - P - P A - 2 -P° w ) - K(P + p A+ 2 -pc w ) 


1 


A U (C - E + P + P a -P ’ 


M(7r 0 7r°) = 

y/2 [ 

= --j= [\ U (T + C — Pi 


ew g Pew) 


yC 


A c{~P — Pa ± Pew ± ~^Pew) 


ew ~ Pew ) — ^ c{Pew ± Pew)] j (17) 


with A u = V ub V* d = AA 3 (p — irj){ 1 — A 2 /2), and A c = V cb V* d = —AA 3 . Throughout the 
present paper, we shall assume the t-quark dominance in penguin type diagrams. In general 
a penguin diagram can be written as 


P — XuPu ± A C P C ± A tPt — —X u Ptu — XcPtc, 


(18) 


where Pt u = Pt ~ P u and Pt c = Pt — P c ■ The t— quark dominance in the Wilson coefficient 
then leads to 3> P c ^> P u and 


P tu ~ Ptr ~P t = P. 


tc 


(19) 


Note that in the presence of large charming penguin |29|, (53J, 154], I55j, l56(, P tc could differ from 


P. This effect can be effective absorbed into inelastic final state interactions(FSIs) and will 
not be discussed in detail here. 
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From the isospin structure of the low energy effective Hamiltonian, the sum T + C and 
Pew + Pew have both isospin 1 = 2. It is then helpful to define 

T = T + C, Pew = Pew + Pew- ( 20 ) 

The ratio R^w is just the ratio of the isospin 1 = 2 part between electroweak penguin and 
tree diagrams. 

In the naive factorization approach bj, these amplitudes have the following typical 
values 


T = 0.9 ~ 1.1, C = -0.33 ~ 0.25, 

P ~ 0 . 1 , P EW = 0.013 - 0.015, 

P% w = -0.0023 ~ 0.003. ( 21 ) 

All the amplitudes are almost real. The ranges of the parameters correspond to the effective 
number of color Nc ranging from 2 to infinity. In the factorization approach, the rescaled 
amplitudes satisfy a hierarchy of 

| / | 3 > \P\ 3 > \Pew\, \PewI ( 22 ) 

which holds for all primed and double-primed amplitudes in nK and KK modes. 

Including the time dependent CP asymmetry, the 7 T 7 t modes provides seven data points. 
A direct fit to the data gives the following best fits corresponding to a local minimum of y 2 
function. 


m = 0.5318;°“ , \C\ = 0.42±g;Sf 1 , 

= -0-84±°;tI , |P|= 0.099±°;°i, 

5p = -0.55 ^; 2 73 7 , 7 = 1-1±S (23) 

with y 2 lin /d.o./ = 0.17/1, where Pew is fixed relative to T through Eq.©. The above result 
show that: 

• The 7T7T data prefer a large \C/T\ = 0.8 ± 0.2, which is in contradiction with the 
factorization based estimation. This is not new, however a large relative strong phase 
of Sc = — 0.84lg;4 7 is also favored. Note that the recent SCET calculation which 
includes charming penguin effects prefers that Im(C/T) should be vanishing at leading 
order j^j]. In the present ttti fit the charming penguin amplitude is not included. The 
considerable uncertainties in the present data us from drawing a robust conclusion 
on the phase of C/T. The situation will be improved when more precise data are 
available in the near future. The large \C\ and its strong phase Sc are required by the 
observed two ratios in Eq. ©. In the following figure (Fig©, the dependences of the 
ratios with \C/T\ and Sc are given. For illustration purpose, we fix other parameters 
to be \T\ = 0.53, |P| = 0.1, Sp = —0.55 and 7 = 1.1, corresponding to their best fits. 

It follows from Fig. [T] that both R + _ and R 00 prefer a large \C/T\ around 0.8. There 
is very little dependence on Sc for i? 0 o- However, the large strong phase of Sc — —1.0 
is required by R+-, namely by the interference between T and C in 7 r _ 7 r° mode. 




FIG. 1: i?oo and R _|_as functions of \C/T\ with different value of <5c.The three curves correspond 

to Sc = —0.5,—1.0,—1.5 respectively. The shadowed bands are the experimentally lcr allowed 
ranges. The other parameters are fixed at their best fitted value in Eq. El 


• The determination of 7 is in agreement with the SM fit. However, in nn system, there 
could be multiple solutions ( 13 , [57|. The x 2 curve as function of 7 is give in FigH 
which also indicates a local minimum of y 2 close to 7 = 0.23. But the corresponding 
best fitted other parameters are \T\ ~ 0.3, | C\ ~ 0.84, 5c — —1.7 and \P\ ~ 0.36 
which looks unreasonable as it favors \C\ \T\ and T ~ |P|. To get rid of the multi¬ 

solutions, one may include the nK modes via flavor SU (3) symmetry. The simplest 
way is to include n + K~ mode only, as it was done in Ref. |2|. The two-fold ambiguity 
in 7 can be easily lifted. 


The value of \P\ ~ 0.1 agrees well with naive factorization estimate in Eq. m while 
T is suppressed. The enhancement of C and suppression of T implies that there could 
be a mixing between a diagram and its color-suppressed counter part. For mr modes, 
it may be due to large FSI through B —> 7 r + 7 r _ ( 7 r° 7 r°) —> 7 r° 7 r°( 7 r + 7 r _ ). A recent 
calculation based on one particle exchange model indeed supports this conjecture 
Such a mixing may also apply to D°n 0 and p°n 0 modes. 
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B. nK modes 

The amplitudes of nK modes are written in a similar way 


9 







A(n + K~) 

A(tT°K°) 

A(n-K°) 

A(tt°K-) 


?c' 


K (T'-P’--P c ew) 


->c 


K(P’ + 3 Pew, 


_ P P' ~ P'ew ~ 3 Pew) ~ K(~P' + P'ew + t^Pew) 

— P + 3 Pew ) - K(P ~ t^Pew)-* 

_ "\/f ^ C P ^ ~ P ~ Pew ~ t^Pew) ~ K(P + Pew + 3 Pew) 


(24) 


with A* = V u bV* s = AA 4 (p — ip), and = V c bV* s = AA 2 (1 — A 2 /2) . Note that in the nK 
modes | A® | is much smaller than |A*|, |A*/A®| ~ 0 . 02 . The suppression of the tree-penguin 
interference and the limited accuracy of the present data make it less effective to extract the 
weak phase 7 from nK modes at the present stage. Thus if one considers 7 iK modes alone, 
it is more useful to take 7 as known from the global SM fit to explore the other hadronic 
decay amplitudes. Taking the SM value of 7 = 1 . 081^21 as biput and also fixing the P' EW 
with the SM relation of Eq. ©, one finds the following solution 


\T'\ = 1.541855 , \C'\ = 2.718:? 1 , 

8c = 3.1 ±0.11 , \P'\ = 0.12 ±0.0023, (25) 

^ _ n 9+0.07 

0p/ — — u.z_ 012 , 


with a x 2 /d-0-f = 2.49/4. From the above result, one arrives at the following observations 

• The 7 iK data favor both large T' and C' with an even larger ratio of |C"/T , | = 1.75±0.7. 
Although the errors are considerably large, it is evident that a large I C'/T'\ ~ 0(1) is 
also favored in nK modes. A similar observation was made in Ref. [5^ where no error 
estimation was given. The large \C'/T'\ is due to the suppression of R n from unity, in 
FigO the value of R n as function of \C'/T'\ is plotted, one sees that in general, a large 
\C'/T'\ with large relative strong phase 8c — 2 can lead to the reduction of the ratio 
Rn- 

• P' is well determined which is about 20 % larger than P, in a good agreement with 
the factorization based estimation that P'/P ~ fx/fn — 1-28. 

• A relatively larger y 2 /d.o./ in nK fit indicates larger inconsistency with the data in 
comparison with that for 7T7 t modes . The sources of inconsistency mainly come from 
Br(7T + K~) and Br(n°K°) and also S(t: 0 Ks)- The best fit prefers a larger value of 
R n ~ 0.9 and a very small S(tt°Ks) — —0.02. 


An alternative way to achieve at a smaller R n is to allow Pew to be larger, which needs 


Ld the SM. Taking 

Pew t 

\T'\ 

_9 7C++12 

— ' '-’-1.38 > 

\c'\ = 

>C' = 

2.76 ±0.28 , 

\p'\ = 

Oi 

7 

11 

- _n f)8 +0 02 
u.uo_ 008 , 

\P'ew 

)pf 
r EW 

— 1 44+ 0 - 08 

— r.^_ 013 . 



1+0.71 

l -0.75> 


(26) 
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FIG. 2: R n as functions of \C'/T'\ with different value of 5c- The three curves correspond to 
5c = 1.0, 2.0,3.0 respectively. The shadowed band is the experimentally la allowed range. 


with x 2 /d.o.f = 0.4/2. Indeed, one sees that a large P EW is favored by the nK data with 
|P/ W /(T / + C ')| = 0.04 ± 0.04. Once P' EW is increasing, the ratio of C'/T' decreases. It 
seems to be a promising way to restore a reasonable value of C'/T'. However, it only holds 
for nK modes. Furthermore, the uncertainties are too large to prevent us to draw any robust 
conclusion on that. The precisions can be improved significantly by making use of the whole 
charmless B decay data connected by flavor SU(3) symmetry. 

It is more difficult to explain the suppression of R in n + K~ and n~K° modes, as both C 
and P E w are absent. If the puzzle of R has to be taken seriously, one needs an enhancement 
of either A or P EW . This possibility will be discussed in sections IIVI and El 


C. KK modes 

For the KK modes, currently only K°K° mode has been observed. The decay amplitude 
is given by 


A(K°K°) = \ z 


- p "+-P‘ 


C" 

EW 


P1)-K{P"--P^PP'a), 


(27) 
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which is dominated by QCD penguin (through b —> d). Neglecting small subleading diagrams 
Pew an d Pa, one can directly extract the amplitude of penguin from the data 

\P"\ = 0 - 2 ±°;f, (28) 


where we have taken the SM value of 7 as input. It is evident that the QCD penguins for 
7 T 7 T, TiK and KK follow a pattern 

1^1 < l-P'l < \P "\, (29) 


and roughly agrees with a factorization based estimation in that the SU(3) breaking effects 
are proportional to the decay constants of the final states, namely 

P' P" Ik 

_ 1^^ _ r^> _ 

P~ P'~ w 

Thus one finds that a separate analysis can indeed provide important information on decay 
amplitudes and test SU(3) symmetry breaking, which can not be obtained by a global y 2 
analysis. 



IV. ANALYSIS USING FLAVOR SU( 3) SYMMETRY 
A. Fit within SM 


We are in the position now to connect the rnr , nK and KK modes through approximate 
flavor SU(3) symmetry. Note that there is no reliable way to estimate the size of the SU(3) 
breaking in theory. From the factorization based approaches the SU(3) breaks in such a 
way that the amplitudes in nK(KK) modes differ from the ones in 71ti{eK) modes by a 
factor of f k / f-K , where fx and f n are decay constants for K and n mesons. There have 
been analysis based on different patterns of SU(3) breaking which could in general lead to 
different results. Here we would like to consider three scenarios of SU (3) relations frequently 
used in the literature: 

Scenario A) All diagrammatic amplitudes for nK^KK) modes are larger than that in 
tht^tiK) modes by a factor fx/C- 


V _C' _P' _ T" _ C" _ P" _ f K 

Scenario B) SU(3) symmetry breaks only in tree diagrams (Til Id 

V T" fx 


C' _P' _ C" _ P" 


Scenario C) Exact SU(3) limit 


rj-if £1/ pf pit pill pH 


(31) 


(32) 


(33) 


As in the first step, the Pew is fixed within the SM through Eq. ©• Thus there are 6 
parameters T, C, Sc, P, Sp and 7 to be fitted by 18 data points. The best-fitted parameters 
as well as the Brs and acps are tabulated in Tab fill 
The results show that: 
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scenario A 

scenario B 

scenario C 

m 

0.52 ± 0.027 

0.51 ± 0.033 

0.52 ±0.032 

\c\ 

0.47 ± 0.042 

0.45 ± 0.053 

0.45 ± 0.053 

6c 

-1.1 ± 0.19 

1+0.21 

—0.19 

1 1+0.21 

0.19 

\P\ 

0.094 ± 0.0014 

0.12 ±0.0019 

0.12 ±0.0018 

6 P 

—0.49±g;?® 9 

—0.45t°;?f 

-0.54;°;^ 

\Pew\ 

0.011 ±0.0011 

0.011 ±0.0011 

0.011 ± 0.0011 

$Pew 

-0.52 ±0.1 

-0.47 ±0.11 

-0.49 ±0.11 

7 

1+0.11 

0.13 

1 1+0.13 

0.17 

1 1+0.14 

U!— 0.17 

xhJd.o.f 

16.2/12 

19.2/12 

21/12 

Br(TT + n~) 

4.7 ± 0.48 

4.8 ±0.62 

4.9 ±0.59 

«cp( 7r+7r_ ) 

0.27 ±0.062 

0.32 ±0.085 

0.37 ±0.096 

Br(ir°Tr°) 

1.7 ±0.31 

1.8 ±0.4 

1.8 ±0.41 

a C p(TT 0 TT 0 ) 

0.36 ±0.11 

0.43 ±0.15 

0.38 ±0.17 

O 

1 

cq 

5.2 ±0.77 

5.1 ± 0.85 

5.1 ±0.86 

a C p( vr vr°) 

0 ± 0.01 

0 ± 0.01 

0 ± 0.01 

Br(ir + K~) 

20 ±0.77 

20 ± 0.84 

20 ± 0.74 

acp(TT + K~) 

-0.1 ±0.02 

-0.097 ± 0.022 

-0.089 ±0.019 

Br(ir°K°) 

9.8 ± 0.49 

9.9 ±0.47 

9.7 ±0.46 

acpi'P 0 K°) 

-0.1 ±0.035 

-0.076 ± 0.03 

-0.068 ± 0.032 

Br(ir~ K°) 

22 ±0.69 

22 ± 0.71 

22 ± 0.68 

Br(ir°K-) 

12 ±0.63 

11 ±0.64 

12 ±0.57 

a C p(n 0 K~) 

0.0055 ± 0.042 

-0.016 ±0.039 

-0.014 ±0.04 

Br(K°K°) 

1.3 ±0.17 

2.3 ±0.35 

0.84 ± 0.13 

Br(K~K°) 

1.3 ±0.17 

2.3 ±0.35 

0.84 ± 0.13 


-0.73 ±0.13 

-0.76 ±0.13 

-0.73 ±0.14 

0 

0 

0.23 ±0.27 

0.51 ± 0.27 

0.52 ± 0.27 

S(n°K s ) 

0.86 ± 0.038 

0.84 ± 0.04 

0.84 ± 0.04 


TABLE II: Global fit to 7T7T, nK and KK modes within SM. The three columns corresponds to 
the three different SU(3) relations used in the fits. 


• The differences among the three scenarios are in general not large. The weak phase 
7 is well determined in all the cases and depends on the SU(3) breaking scheme 
very weakly. All the three scenarios give 7 ~ 1.1 in a good agreement with the SM 
value with differences less than 10 %, which manifests that 7 can be reliably extracted 
using the diagrammatic approach. The Xmin curves as function of 7 are given in 
Fig 0, Comparing with the one from tttt fit, one finds a significant improvement on 
the precision of 7 determination. The three patterns lead to roughly the same \T\ and 
\C\ with \C/T\ ~ 0.8. Note that for small \C\ we find no consistent fit. For example, 
if \C\ is fixed at C — 0 . 2 , a very big Xmin/d-o.f = 44.6/12 is obtained. The major 
difference is the value of \P\. The scenario B) and C) prefer a \P\ which is ~ 20% 
larger. 


13 







• Among the three cases, the scenario A) that all the primed (double-primed) amplitudes 
are larger than the unprimed (primed ) ones by a factor of fx/fn gains the lowest 
X 2 /d.o.f = 16.2/12, which indicates a better consistency in comparison with the other 
two. The exact SU(3) scenario gains the largest y 2 /d.o./ = 21/12, which clearly 
indicates that the flavor SU (3) symmetry in charmless B decays must be a broken 
one. 

• The main source of the inconsistency comes from the Br(TT + K~), Br^^K 0 ) and 
Br(n~K ). The best fits in scenario A prefer a larger Br(n + K~) ~ 20, a small 
Br(n°K ) ~ 9.8 and a small Br(n~K ) ~ 22. Namely, within the current parame¬ 
terization, it is not possible to arrive at the observed ratios R n and R. Thus the ttK 
puzzles remain. 

• For the predictions for the KK modes, the scenario A gives Br(K°K °) = 

Br(K~K°) = 1.2, while the other two give 1.7(scenario B) and 0.84 (scenario C) re¬ 
spectively. The branching ratio of K + K~ is predicted to be zero and all the predicted 
direct CP violation are vanishing, due to the lack of interferences between amplitudes. 

Other possibilities of SU(3) breaking include the SU(3) breaking in strong phases, which has 
been discussed in Refs. |34j . The current data favor a small SU(3 ) breaking in the strong 
phase of QCD penguin. This breaking effect can significantly modify the correlation of direct 
CP asymmetries between 7T7 t and tiK modes Q. 


B. Fit with free electroweak penguin Pew 


In the next step, we consider the possibility that Pew is free from the SM constraint. 
The fitting results with free Pew under three scenarios are given in Tab lllll 

In this case, one finds that all the fits prefer a larger value of \Pew\ = 0.23 ~ 0.29 with 
a large strong phase Spew ~ 0- 6 ~ 0.7 relative to T as the corresponding best fit of T has 
a strong phase of —0.5 ~ —0.4. A large Pew with a large strong phase relative to T can 
naturally explain the suppression of R n and also R 2 [1, 3 EM 20l E, 37, -38. 39]. Naively 
speaking, all ttK modes are QCD penguin dominant, the ratios R„ and IP should be very 
close to unity. The corrections arise from either tree type diagrams or electroweak penguins. 
The former is CKM suppressed in eK modes and is constrained by nn data. Thus an 
enhancement of electroweak penguin is needed. 

As the two ratios R n and R 2 are both P E w sensitive, they can be used as probes of 
electro-weak penguins. We parameterize the deviation of SM by introducing a complex 
parameter k 


Pew 

f 


E> 

It 


.SM 

EW 


K. 


(34) 


In FigEl the two ratios are plotted with different values of |k| and its strong phase 5k. For 
demonstration reasons, the other parameters are fixed at the best fitted values in SM in the 
first column of Tab. |TT] according to the scenario A, namely 


T = 0.52, C = 0.47, 5 C = -1.1, P = 0.094 and 5 P = -0.49. (35) 
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FIG. 3: Xmin as functions of weak phase 7. The three upper curves correspond to the three fit in 
TabUll The solid, dashed and dotted curves correspond to scenario A, B and C respectively. The 
lower curve (dot-dashed) corresponds to the fit only to 7T7 t modes in Eq. (fTTl) . The shadowed band 
is the allowed range from global SM fit. 


The figures indicate strong dependences on 5k and |ac| for both ratios. In the SM case, i.e 
|ac|=1 and 5 k = 0, R n is expected to be above 1.0 in contrast with the experiment. For |k| = 1 
and large S K ~ 1.5, R n is reduced but R2 is enhanced and departs away from the allowed 
range of the data. Thus to simultaneously explain both measurements, one needs a large 
\k\ ~ 2.0 ~ 3.0 and a large strong phase of 5 K ~ 1.0 ~ 1.5. This observation is confirmed by 
the global fits with Pew f ree in Tab. 11111 which gives 


Pew 

f 


(3.1 ± i.3) e *( 1 '0 2 ±o.5) x -^g-2 (scenario A) 
(4.8 ± i. 5 ) e *( 1 - 06 ±o. 53 ) x -^g-2 (scenario B) 
(3.1 ± i. 5 ) e *( 1 - 03 ±o. 5 ) x iq-2 (scenario C). 


(36) 


With Pew being free , the best fitted n°K° mode is found to be Br(ir°K 0 ) — 11 ± 1.9 in 
scenario A, in a remarkable agreement with the data. The central values of the two ratios 
are found to be R n ~ R 2 — 0.9. 

The enhanced Pew with large strong phase will result in different predictions for the CP 
asymmetries. In most decay modes the predicted acps are much smaller |2], However, the 
most important one is Ocp( 7 r _ 7r°) which should be exactly zero in SM. But now it prefers 
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scenario A 

scenario B 

scenario C 

m 

0.52 ± 0.028 

n ^1+0.037 
u ' oi -0.045 

0.51 ±0.035 

\c\ 

0.45 ± 0.052 

0 44+0.096 
u '^-0.062 

0 44+0-078 
u ’^-0.064 

Sc 

—0 Qfi+ 0 ' 23 

U.JO_0.21 

-0.92 ±0.25 

-0.98 ±0.24 

\P\ 

0.093 ±0.0015 

0.12 ±0.0019 

0.12 ±0.0019 

Sp 

—0 52 +01 
u ' oz -0.13 

— 0 49+ 0 - 1 

U ' 4ly -0.24 

-o. 59 + 0:22 

1 Pew\ 

0 noQ+0-0096 
u.uzo_ 0 on 

0.027 ± 0.014 

0.029 ±0.013 

Spew 

U.OU„o.4l 

p\ '7+0.23 

U ’ * —0.35 

0 (54+ 0 - 23 
U.OO_o.32 

7 

1+0.13 

—0.18 

i 1+0.16 

0.36 

1+0.17 

0.26 

X 2 miJ d -0-f 

13.2/10 

15.9/10 

18/10 

Br(TT + n~) 

4.7 ±0.53 

4.8 ± 0.84 

4.9 ±0.7 

acp(p + P~) 

0.29 ±0.077 

0.34 ± 0.14 

0.39 ±0.14 

Br(ir°Tr°) 

1.6 ±0.38 

1.7 ±0.66 

1.7 ±0.58 

a C p( 7T 0 7T 0 ) 

0.14 ±0.18 

0.16 ±0.26 

0.12 ±0.26 

1 

5.4 ±0.89 

5.3 ± 1.2 

5.2 ±1.1 

acp (vr 7T°) 

-0.085 ± 0.045 

-0.11 ± 0.061 

-0.11 ±0.06 

Br(n + K ~) 

20 ±0.85 

20 ± 1.1 

20 ± 0.84 

acp{n + K~) 

-0.11 ±0.024 

-0.1 ±0.035 

-0.093 ± 0.027 

Br(TT°K °) 

11 ± 1.9 

11 ± 1.8 

11 ± 1.9 

a C p(TT°K°) 

-0.034 ± 0.045 

-0.027 ± 0.043 

-0.02 ± 0.042 

Br(ir~ K°) 

22 ± 0.73 

22 ± 0.74 

22 ± 0.73 

Br(TT°K~) 

12 ±2.2 

12 ±2.2 

12 ±2.2 

acp(p°K~) 

0.033 ± 0.059 

0.012 ± 0.063 

0.01 ±0.057 

Br(K°K°) 

1.3 ±0.21 

2.3 ± 0.56 

0.81 ±0.17 

Br(K~K°) 

1.3 ±0.21 

2.3 ± 0.56 

0.81 ±0.17 

) 

-0.7 ±0.17 

-0.73 ± 0.22 

-0.71 ±0.2 

O 

0.4 ±0.31 

0.54 ± 0.41 

0.55 ±0.36 

S(ir°K s ) 

0.86 ±0.039 

0.84 ± 0.043 

0.84 ± 0.042 


TABLE III: The same as TabUll but Pew is taken as a free parameter to be determined directly 
from the data. 

a negative value of acp( vr _ 7r°) ~ —0.1, which is in agreement with the current preliminary 
data of acp = —0.02 ± 0.07 and can be examined with more precise data in the near future. 

It needs to be emphasized that the large Pew here means a relative enhancement to tree 
type diagram T, not to QCD penguin diagram. It was claimed recently in Ref. Q that 
there was no clear indication of large Pew /-P, which does not necessarily contradict with 
the conclusions in the present paper. Furthermore, using the numerical value of T and C 
obtained in Ref. 0. we find a similar result as in Eq. ©■ Note that the ratio Pew/P is 
subjected to large theoretical uncertainties, it is better to use the values relative to T for 
exploring the electro-weak penguin and new physics as it is free from hadronic uncertainties. 
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FIG. 4: Xmin as functions of weak phase 7. The three upper curves correspond to the three fit in 
Tab II I II The solid, dashed and dotted curves correspond to scenario A, B and C respectively. The 
shadowed band is the allowed range from global SM fit. 


C. Effects of color suppressed electroweak penguin P EW 

In the previous discussions, the color-suppressed electroweak penguin diagram P EW 
neglected. However, among all the subleading diagrams P EW * s the on ly one giving con¬ 
tribution to high isospin state / = 2(3/2) in Tnr(nK) modes. Furthermore, it cancels the 
isospin / = 0(1/2) components of Pew just like C cancels that of T and directly contribute 
to the ratio in Eq. ©■ Since the current data indicate a sizable C, it is still possible that 
there will be an enhancement of P EW as well. In the SM, Pew is fixed relative to T. How¬ 
ever, given a large relative strong phase i.e negative interference between P EW and P E w 
a large value of P EW is possible within the SM. Note that the 7 t°K° mode depends on 
Pew + P E w = Pew/3 + 2P E w/3- Even the first term is constrained by Eq. ©, the second 
term can still enhance the decay rate of 7 t°K° mode without violating the SM relation. The 
similar arguments also applies to ti°K~ modes which depends on P E w + 2 P E w /3- 

To see the effects of P E wi we give iu FigEl the ratios of R n and R 2 as function of P EW and 
its strong phase <5pc^ ■ In the numerical calculations we take the values of other amplitudes 
from the SM fit, according to the first column of TabUTlor Eq. f|35[) . The values of P EW and 
its strong phase are automatically generated from Eq. ©• It follows from the figure that for 
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FIG. 5: R n and R 2 as function of 5 K with different value of |k|. The three curves correspond 
to |k| = l.O(solid), 2.0(dashed, 3.0 (dotted). The shadowed bands are the experimentally allowed 
ranges. Other parameter are fixed at the SM best fitted values (column A in Tab DU). 


a small | Pp jW | ~ 0.1 the predicted value of R n is far above the data for all values of 3 P c w - 
To account for the current data | Pp W has to be greater than ~ 0.4. The strong phase 8 P c 

receives strong constraint from R 2 . For | Pp W | in the range 0.1 ~ 0.6, a large negative value 
of 5 P c = —0.2 ~ —2.5 is favored. 

r EW 

Including Pp W as a new free parameter while keep Eq. S. we find 

I did ~ 0.025 ± 0.021, S P c w = -2.24«;1 (37) 


and 


Pew 


0.03 


+0.02 

-0.013) 


W — 0-52^0^2) 


(38) 


with a x 2 fd.o.f = 11.3/10. The best fits of other parameters are listed in the first column 
of Tab IIVI Clearly, there is a strong cancellation in the sum of Ppw + Ppw as required by 
Eq.©. It is of interest to note that for both tree and penguin diagrams the color suppressed 
diagrams are not necessarily suppressed. Furthermore, the current data suggest that 


c 


P c 

r EW 

T 


Pew 


(39) 


Thus the relative enhancements are likely to be universal. This is again in favor of the 
conjecture that it has a strong interaction origin which is flavor independent. Comparing 
with Eq. (gup , one sees that Pp W is on the lower side to account for the suppression of R n . 
The best fitted ratios are R„ ~ 0.89 and R 2 — 0.78 respectively. Thus a large Pp W improves 
the agreement with the data. 
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FIG. 6: Ratios of R and R 2 as functions of <5pc and \P EW \- The three curves correspond to 
\Pp jW \={).{)\ (solide), 0.04 (dashed) and 0.06 (dotted) respectively. Other parameters are fixed at 
their best fitted value in SM (according to column A of Tab IH. The shadowed bands are the 
allowed range by the current data. 


Taking both P EW and P EW as independent free parameters, we get the following fit result 


P£ w \ = 0.016 ± 0 . 02 , 


5 P c = -2.59+°^ 
r — 1 ./ 


(40) 


and 


IP) 


EW\ 


— 0 027 +0 ' 016 

— u.uz/_ a014 , 


S Pbw = 0.69±°;1, 


(41) 


with Xmin/d- 0 -f = 6.6/8. The color suppressed electroweak penguin is found to be reduced 
but still significant. In this case, all the nK ratios are well reproduced. Note that the 
best fits correspond to \Pew/T\ = 0.032 ± 0.018 which again implies a violation of the SM 
relation. 

Finally, we emphasize that the large P E w within the SM may distinguish itself from the 
one beyond the SM by the prediction of direct CP violation in B —> 7r _ 7r° which should be 
vanishing in the former case and small but nonzero in the later. 

Given a small P EW relative to P E w, the current data may imply new physics beyond 
the SM. New physics models significantly contributing to charmless B decays may include 
various SUSY models |2^|, Z' models from extra U( 1) gauge symmetry |21, 0] and two- 
Higgs-doublet models (2HDMs)(60). Among various versions of two-Higgs-doublet models, 
the general 2HDM with spontaneous CP violation can provide rich sources of CP violation 


6lL 1621.163L |6J, |65( and significant corrections to the electro-weak penguin through charge 
or neutral Higgs boson exchanges. The models with 4th generation may also give sizable 
contributions. For model with both two-Higgs-doublet and 4th generation quark, the effects 
could be more significant through neutral-Higgs loops0, [(37, [68). 
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V. NONFACTORIZABLE DIAGRAMS 


We now go a step further to discuss the effects of other subleading nonfactorizable dia¬ 
gram such as E, A and Pa- They are expected to be very small from factorization based 
estimations. However, in the presence of large FSI, there could be mixing among diagrams 
which may enhance the sizes of subleading diagrams fpi Efi 10 . In view of the current 
puzzling pattern of the data, the possibility of anomalously large nonfactorizable diagrams 
can not be excluded j|69j]. 

Due to the limited number of data points, it is not possible to directly extract all of them 
simultaneously from the current data of 7T7t, ttK and KK . Instead, to obtain the typical 
sizes of those diagrams, we shall consider several typical cases, in each case only one of the 
diagrams is assumed to be dominant while the other two are small. For simplicity, we only 
consider the case where Pew is fixed as the SM value. 


A. IF-exchange diagram E 


It has been argued that sizable E with constructive (destructive) interference with C(T) 
can help to understand the large value of \C/T\ obtained when E is absent 0,0,0. Since 
the main contribution to 7r + 7r - is from T + E while it is C — E for 7r°7r° mode. To illustrate 
how E improves the consistency with the 7T7 t data with a small \C/T\, we fix the tree and 
color-suppressed tree to be 


\T\ = 0.9, \C\ = 0.3, and S c = 0, (42) 

respectively. The QCD penguin is fixed at \P\ = 0.1, Sp = —0.5 and Pew is fixed at its 
SM value. In this case, the dependence of R 00 and i? + _ with E and its strong phase Se is 
plotted in Fig [7] 

As shown in the figure, for a small ratio of |C/T|=0.3, the data require a large \E\ ~ 
0.3 ~ 0.5 with a large strong phase of ~ 2.0. However, The origin of large E is still a 
challenge for theory. 

Assuming E is dominant, we find a big value from a fit to the data 

\E\ = 0.461°;“ , S E = 2.86l°;£. (43) 

The whole fit result is given in the second column (B) of Tab If VI The Xmin/d-°-f is 13.2/10. 
Note that although \C/T\ is reduced to ~ 0.28, the relative strong phase Sc is found to 
be large Sc = — 2.l!2’72 in contradiction with factorization based estimates. Thus only 
introducing a large VF-exchange diagram will not be enough to coincide with factorization 
results. 

An alternative way to reduce C/T is to make Pt u significantly different than P or P tc . 
This is less likely as P u is greatly suppressed relative to P t by small w-quark mass in Wilson 
coefficient. The important FSI process such as charming penguin only affects P tc . A fit 
taking P tu and P tc as independent parameters shows that the best fitted P tc is close to the 
factorization estimate of P while P tu is large and compatible with T in size 0], which is 
quite unreasonable. 

As it was mentioned previously, an enhanced E has no effect in ttK modes, thus can not 
solve the nK puzzle and explain the obtained large C'/T'. The best fitted branching ratios 
in nK still exhibit the puzzling patterns. 
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FIG. 7: i?oo and R-\ as functions of E with different values of 5e The three curves correspond to 

5e =1.8 (solid), 2.0 (dashed), 2.2( dotted). The shadowed bands are the experimentally lcr allowed 
ranges. 


B. Annihilation diagram A 

The annihilation diagram A appears in n~K° and tt°K~ modes. The current data indicate 
that both of them are large in comparison with tt + K~ mode, which is characterized by the 
suppression of R 2 and R. This may imply a sizable A in these modes as discussed in the 
previous sections. In Fig|H] the two ratios are given as functions of A and 5 a- The other 
parameters are fixed at the SM best fit value in Tab. ini- 

One sees from the figure that the two ratios have a similar behavior under the changing 
of A and 5a- For A ranges between 0.3 ~ 0.8, both ratios decrease with 5a growing. For 
\A\ = 0.5 ~ 0.8 the curves fall down to the experimentally allowed ranges at 5a > 2.0. Since 
both ratios contain the same term A + V, they have similar dependences on A and 5a- This 
will lead to a cancellation for the ratio of the ratios: R/R 2 =2Br(n 0 K~) / K°) should 
be very close to unity, which agrees with the data well. 

Assuming A is dominant over other subleading diagrams and using the scenario A for 
SU(3) breaking, we find 


M = 0.231" H , S A = 2.77l»S, (44) 

with Xmin/d-O-f = 13.8/10. The whole fitting results are listed in column (D) of Tab IIVI 
One sees that the best fitted value of A is moderate, which helps to reduce R and R 2 but is 
not large enough to reproduce the central values of the two ratios. Note that even this value 
of |A| ~ 0.2 is still much larger than that from factorization estimation which is suppressed 
by a factor of fs/^B- 
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FIG. 8: R and R 2 as function of 5a with different values of A. The three curves correspond 
to A =0.3 (solid), 0.5 (dashed), 0.8(dotted).The shadowed bands are the experimentally allowed 
ranges. 


C. Penguin annihilation diagram Pa 


In 7T7T modes, the penguin induced annihilation diagram Pa contributes to only low isospin 
final states (/ = 0) and it always comes together with QCD penguin P in nn modes. 
Although Pa is often neglected in the literature, its effects are however effectively absorbed 
into QCD penguin in nn modes. Thus the QCD penguin extracted from mr modes are 
effectively P = P + Pa- In general Pa acquires a strong phase different than that of P, 
namely the strong phase of P and P are different. Since there is no SU(3) counter part 
P' A appearing in 7 iK modes, P' is still have the same strong phase as that of P in SU(3) 
symmetry. This introduces an effective SU(3) breaking in strong phase between P extracted 
from mi modes and P' from nK. A fit to the current data gives the following values 

P A = 0.035lo!oi5 , bp A = -2.26 ± 0-48. (45) 

Thus its size is compatible with that of electro-weak penguin. The Xmi n /d-o.f is found to be 
14.3/10. The best fitted other parameters can be found in the column (C) of Tab IIVI The 
best fitted value corresponds to P ~ 0.093 and dp ~ —0.85, which corresponds to a strong 
phase shift of 


A 5p = 5'p — Sp — 0.38. 


This is in a good agreement with the previous analyses in 


Refs. 0,0. 


(46) 


VI. IMPLICATIONS FOR KK MODES 

The KK modes are much more sensitive to the subleading diagrams 8, A and Va- The 
K°K° mode is dominated by QCD penguin as discussed in section ITTT1 The decay amplitudes 
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A 

B 

C 

D 

\T\ 

0.51 ±0.027 

n q4 +°- 31 

U - J ^~0.24 

0.53 ± 0.026 

0.51 ± 0.027 

\c\ 

0.47 ± 0.048 

0 26 +0 ' 2 
u - ZD -o.n 

0-^iao46 

0.48 ± 0.044 

Sc 

—0 Q7 + 0 - 28 

U.y/_ 0.22 

9 -1 -(-0.83 
z * 1 -0.72 

— 1 i +°- 2 

— 1 1+°' 2 

\P\ 

0 0Q4 +U ' uui8 
u.uau_ 0002 | 

0.097 ± 0.0022 

0.094 ± 0.0015 

0.093 ± 0.0015 

6 P 

-0.41 ± 0.14 

—0 2 fi +0 ' 087 

u ' zd -0.13 

-0.47+^f 

— 0 c >4 +0 - 1 
U ' O4 -0.13 

1 Pew\ 

0 03+ 0 - 019 

u - uo ~0.013 

0.011 ± 0.0011 

0.011 ± 0.0011 

0.011 ± 0.0011 

Sp EW 

0-52±g;fi 

—0 28 +0 ' 14 
u ' zo -0.16 

—0 54 +0,11 

u ' jy: -0.25 

-0.53 ±0.11 

1 pc 1 

\ r EW\ 

o.o25;°;°2 8 

0 (fix) 

0 (fix) 

0 (fix) 

SpC 

r EW 

—2 24 +0,21 
z - zq: -0.61 

0 (fix) 

0 (fix) 

0 (fix) 

\E\ 

0 (fix) 

0-46±g;if 

0 (fix) 

0 (fix) 

5e 

0 (fix) 

2-86±8;17 

0 (fix) 

0 (fix) 

\Pa\ 

0 (fix) 

0 (fix) 

0.035+^f 

0 (fix) 

Sp A 

0 (fix) 

0 (fix) 

-2.26 ±0.48 

0 (fix) 

14 

0 (fix) 

0 (fix) 

0 (fix) 

0 94 + 0 - 12 
U - ZO -0.087 

5a 

0 (fix) 

0 (fix) 

0 (fix) 

9 77 '+ 0.28 

z -' ' —0.35 

7 

n qq+0-13 
U - JO ~0.16 

n 04 + 0.12 

n oq+o.13 

U.JO —0 43 

n oq+ 0-12 
u - yj -0.15 

X-min/ d'O' f 

11.3/10 

13.2/10 

14.3/10 

13.8/10 


TABLE IV: Global fits to hit, nK and KK data with subleading diagrams. The four columns 
corresponds to the four cases in each of them one subleading diagram is set to be free parameters 


of the other two modes read 

MK'K i = - - I'Tl - A../'"], 

A(K-k<‘) = A„(-P" + 1 PS' W + A") - A C (P" - \p% w ). (47) 

The K + K~ modes depend only on the subleading diagrams E and Pa- Thus it provides an 
ideal avenue to explore their effects. From the current upper bound of Br(K + K~) < 1.8 |28j|. 
and the SU(3) relation of scenario A, the exchanging diagram E receives a constraint of 

\E\ <0.3, (48) 

which limits its contribution to 7T7 t modes to be moderate as the best fit to mr modes require 
an \E\ ~ 0.48 in Ea. (14311 . It is expected that stronger constraint on E will be found with 
more precise data in the near future. 

The direct CP violations for K + K~~ and K~K° reads 

a C p(K+K-) 

21 X u X c 11 P^E" | sin 7 sin 5 

“ \\ U \ 2 (\E"\ 2 ± \P'{\ 2 - 2 \P'{E"\ cos5) + \\ C \ 2 \P'{\ 2 + 2|A U A C ||P/(| cos 7 (|P"| cos5 - \P%\) 

(49) 
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and 


a CP (K-K °) 


" |A tt | 2 (|A"| 2 + |P"| 2 


2\\ U \ C \\P"A"\ sin 7 sin 5' 

2\A"P"\ cos5') + |A c | 2 |P "| 2 + 2 |A u A c ||P /, |cos 7 (|A ,/ |cos 5 / - |P"|) 

(50) 


where 5 = 5 p ^—5e" and 5' = 5 P " — 5a"- In the expression of acp(K~K° ) the color suppressed 
electroweak penguin are neglected. A nonzero acp(K + K ~) will definitely indicate both 
nonzero E and Pa- In spite of the small branching ratio, in the case that E and Pa are 
compatible in size, and the strong phase difference is large, then the direct CP violation 
could be significant. 

Unlike in the nK modes where A is suppressed by a factor |A®/A®| = 0( A 2 ), in K~K° 
mode, it is not suppressed. Thus it is promising to probe A in K~K° mode. A sizable 
annihilation diagram A will show up either through the difference between Br(K°K° ) and 
Br{K~K°) or through the nonzero direct CP violation i.e. acp{K~K°) ^ 0. 




FIG. 9: Branching ratio and direct CP violation of K~K° as functions of A with different values 
of 5a- The three curves correspond to 5a =0 (solid), 2.5 (dashed), 3.0(dotted) respectively. The 
shadowed band is the experimentally la allowed range. Other parameters are fixed at their best 
fitted value in SM (according to column A of Tab HU). The SU(3) breaking effects are taken into 
account according to scenario A. 


In FigEl the decay rate and direct CP violation of K~K° mode are plotted as function of 
A. In the numerical calculations, the value of A” and P" is calculated from the best fitted 
value of P according to scenario A. For 5a in the range 0 ~ 2.5, both decay rate and direct 
CP violation increase with A increasing. One sees from the figure that for |A| = 0.15 ~ 0.2 
and 5a = 2.5 ~ 3.0, acp{K~K°) can reach 0.2 ~ 0.4 with the branching ratio in agreement 
with the current data. Thus with significantly large subleading diagrams, it is promising to 
observe large direct CP violation in this decay mode. 
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VII. SUMMARY 


In summary, we have systematically studied charmless B decays B —> tttt , ttK and KK, 
following a strategy that making independent analysis for 7 r 7 r, ttK, KK modes individually 
as the first step and then connecting them through various SU(3) relations. The separated 
analysis allowed us to clarify the origins of the inconsistence or puzzles revealed by the 
current data. Independent analysis on nn and nK modes both favor large ratio of C/T and 
C'/T' with large strong phases, which suggests that they are more likely to originate from 
long distance strong interactions rather than large non-factorizable exchange diagrams E. 
The sizes of QCD penguin diagrams in 7 T 7 t, itI\ and KK are independently extracted and 
were found to follow a pattern of SU(3) breaking in a good agreement with factorization 
estimation P'/P ~ P"/P' ~ fx/fn- Global fits to these modes have been carried out 
under various scenarios of SU(3) relations. All the results show good determinations of 
weak phase 7 in consistency with the Standard Model (SM) and prefer a large clectroweak 
penguin {Pew) relative to T + C with a large strong phase. Within the SM, it may require an 
enhancement of color suppressed electro-weak penguin {P^w) with destructive interference 
to Pew ■ The possibility of the presence of new physics effects can not be excluded. We 
have also investigated the possibility of sizable contributions from nonfactorizable diagrams 
such as E, A and Pa- Their sizes could be significantly larger than the expected ones. The 
typical sizes of \E\ could reach to \E\ ~ 0.3 as required by the tut and KK data, |A| could 
reach to \A\ ps 0.2 while \Pa\ has a typical value of \Pa\ ~ 0.03. The sizable subleading 
diagrams may change significantly the predictions for the yet to be seen K + K~ and K~K° 
modes. The CP violation in K~K° modes could reach a,cp — 0.2 ~ 0.4 for a large value of 
A and 5a in the range of |A| = 0.15 ~ 0.2 and 5a = 2.5 ~ 2.5. It would be encouraging 
to expand the investigation of subleading diagrams to decay modes involving r) and 7 / final 
states. Although these decay modes receive significant contributions from additional flavor 
singlet penguin diagrams pol. Ful Fil. j73| or nonstandard contributions through cc [?4|. [a| 
or QCD anomaly j76|, [77[, as more data points are involved, stronger constraints on the 
subleading diagrams will be expected. Thus it will enable us to test the SM using the full 
diagrammatic decomposition. 
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